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ABSTRACT 

Dust grains can be efficiently accelerated and shattered in warm ionized medium (WIM) 
because of the turbulent motion. This effect is enhanced in starburst galaxies, where gas is 
ionized and turbulence is sustained by massive stars. Moreover, dust production by Type II 
supernovae (SNe II) can be efficient in starburst galaxies. In this paper, we examine the effect 
of shattering in WIM on the dust grains produced by SNe II. We find that although the grains 
ejected from SNe II are expected to be biased to large sizes (a > 0.1 |xm, where a is the grain 
radius) because of the shock destruction in supernova remnants, the shattering in WIM is ef- 
ficient enough in ~ 5 Myr to produce small grains if the metallicity is nearly solar or more. 
The production of small grains by shattering steepens the extinction curve. Thus, steepening 
of extinction curves by shattering should always be taken into account for the system where 
the metallicity is solar and the starburst age is typically larger than 5 Myr. These conditions 
may be satisfied not only in nearby starbursts but also in high redshift (z > 5) quasars. 

Key words: dust, extinction — galaxies: evolution — galaxies: starburst — H II regions — 
supernovae: general — turbulence 



1 INTRODUCTION 

Type II supernovae (SNe II) are considered to be one 
of the grain production s ources in the Univers e (e.g . 
Kozasa. Hasegawa, & Nomotd 1 19891 : iTodini & Ferraral 1200 ll : 
Nozawa et al.l 120031) . The significance of SNe II in the grain 



production is enhanced if the cosmic age is so young (typically 
at redshift z > 5) that low-mass stars, i.e. asymptotic giant 
branch (AGB) stars and Type la supernovae, cannot contribute 
signific antly to the dust form ation dDw ek. Gallia no! & Jonesl2007l : 
but see IValianteetaT] 12009). or if the current starburst is strong 
enough. For some nearby blue compact dwarf galaxies (BCDs) , 
the latter condition may be satisfied dHirashita & Hunt! 1 20041 : 
Take uchi et alj2005h . Thus, the dust production by SNe II is tested 
by high z objects and nearby starbursts. 

In galaxies where active star formation (starburst) is occur- 
ring, dust grains are not only produced and ejected from stars 
but also processed in the interstellar medium (ISM). In such star- 
forming galaxies, it is expected that the supernova (SN) rate is 
enhanced, which leads to an efficient destruction of large grains 
with a > 0.1 urn, where a is the grain radius, by SN shocks 
(e.g jjones et alJl994l ; |jones. Tielens. & Hollenbachlll996l) . Also a 
large amount of ionizing photons are supplied and H II regions de- 



velop. Indeed gi ant H II regions w ith sizes of > 100 pc are found in 
nearby galaxies dKennicuttl l 1 984 ). Such lar ge H II regions can als o 
be expected theoretically with a starburst dHirashita & Huntll2004h . 
Moreover, the size is larger than the typical expansion radius of SN 
shell s when the dust cond ensed in SNe II is finally supplied to the 
ISM dNozawa et"al]|2007l . hereafter N07). Thus, it is reasonable to 
consider that the dust ejected from SNe II is supplied to the ionized 
regio ns in starburst galaxies . 

iHirashita & Yanl d2009l hereafter HY09) show that shattering 
occurs efficiently in warm ionized medium (WIM), where grains 
are ef ficiently accelerated by magnetohydrodynamic (M HD) turbu- 
lence ( lYan & Lazarianll2003l : lYan. Lazarian. & Draindl2004r) . The 
same mechanism is also expected to work in the ionized regions 
of starburst galaxies. Indeed, turbulence is ubiquitous in the ISM, 
and the collective effects of OB stellar winds and supernovae 
(SNe) can play an impo rtant role in sustaining turbulence (e.g. 
lElmegreen & Scaldl2 004 ) . Therefore, it is probable that the grains 
ejected from SNe II into the ionized regions are efficiently shattered 
by turbulence. 

Since large grains with a > 0.1 \im are suggested to be in- 
jected into the ISM from SNe II selectively due to the destruc- 
tion in hot plasma produced by the reverse and forward shocks in 
supernova remnants (SNRs) (N07fl shattering in WIM is impor- 
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We call this destruction process 'destruction in SNRs'. 



2 Hirashita et al. 



tant for the production of small grains. If a significant amount of 
small grains are produced, optical-ultraviolet (UV) grain opacity 
is enhanced and the slope of the extinction curve becomes steep. 
In particular, the extinction curves of young starbursts are used to 
constra in the composition and size distribution of grains formed in 
SNe II dMaiolino et al.l2004bl : lHirashita et al.l2005t> . Thus, it is im- 
portant to quantify the effect of grain shattering on the extinction 
curve. 

Also in the observational context, production of small grains 
in starburst environment is worth investigating. The modification of 
grain size distribution should have an impact on t he dust extinction 
and emission properties (e.g. iDopita et alj|2005h . By studying the 
spectral energy distributions (SED s) of dust and stars o f some ac- 
tively star-forming dwarf galaxies. lGalliano et al. (2005) show that 
the grain size distribution is biased to small sizes (~ a few nm). The 
extinction curves of starburst galaxi es in general sh ow a significant 
reddening in the optical-UV range ([Calzetti 200 1]), indicating that 
there s hould be some contribution from small grains. iGalliano et aU 
(2005) also suggest that their results are consi stent with shatter- 
ing a nd erosio n of ISM grains by SN sh ocks d Jones et ail 1 19941 
1 19961 : see also iBorkowski & Dwek|[l995r) . Both shock and turbu- 
lence are efficient drivers of the relative motion between grains, 
since the grain acceleration occurs in a way strongly dependent on 
the grain size JShulllll977l : llVlcKee et alj [l987; Jone s et alj[l996l ; 
lYan et al.l2004l) 3 However, at present, it is not known which of the 
two drivers are more important. Thus, in this paper, as a first nec- 
essary step to fully understand the possible mechanisms of small 
grain production, we focus on turbulence as a possible driver of 
grain shattering. 

This paper is organized as follows. We explain the method 
used in this paper in Section [2] and describe some basic results of 
our calculations in Section [3] We discuss the results and mention 
some observational implications in Section [4] Finally, Section [5] 
gives the conclusion. 



2 METHOD 

We consider a young starburst, where dust is predominantly sup- 
plied by SNe II. We then calculate the modification of grain size 
distribution by shattering in WIM. Finally, we examine if the shat- 
tering effect is apparent in the extinction curve or not. In this sec- 
tion, we first explain the initial grain size distribution in the ejection 
from SNe II (Sections l2.1l and l2~2t . Next we rev i ew the treatm ent of 
shattering which was adopted from lJones et alj l|1994, 1996) in the 
previous paper (Section [2~3l >. We use the grain velocities calculated 
by a MHD turbulence model to obtain the relative grain velocities 
in shattering (Section [2T4l >. We also describe the calculation method 
of extinction curve (Section|23](. Throughout this paper, grains are 
assumed to be spherical with radius a. 



2.1 Initial grain size distribution 

The size distribution of grains ejected from SNe II into the ISM 
(WIM) is adopted from N07. This size distribution is used as the 
initial condition for the calculation of shattering in WIM. N07 
have shown that the size distribution of dust formed in the ejecta 
is largely modified by the destruction by sputtering in ionized gas 



2 Galactic-scale bulk motions such as collective outflow from stellar feed- 
backs are not efficient in producing the relative motions among grains. 



heated by the reverse and fo rward shocks. Not only N07 but also 
Bianchi & SchneideJ d2007l) treat the effect of shock destruction 
in SN e II. N07 consider some aspects which iBianchi & Sch neider 
d2007l) did not take into account: N07 solve the motion of dust 
grains by taking into account the gas drag and treat the destruc- 
tion of dust in the radiative phase as well as in the non-radiative 
phase of SNRs. Thus, we adopt N07's results for the size distribu- 
tion of grains ejected from SNe II into the ISM, although qualitative 
behaviour o f our results are common ev en if we adopt the size dis- 
tribution of Bianchi & Schneider (2007). Also we should mention 



Bianchi & SchneideJ d2007l) neglect the effects 



that both N07 and 
of dust electrical charge and the effects of magnetic fields on grain 
kinematics. These physical processes, which should be quantified 
in future work, are further discussed in Section l4~4l 

Here we briefly summarize the calculati on of N07. N07 starte d 
from the grain size distribution calculated bv lNozawa eTall d2003h . 
who treat ed dust nucleation and g rowth in SNe II based on the SN 
model of lUmeda & Nomotol [ 2002). Then, N07 took into account 
the dust destruction by kinetic and thermal sputtering in hot gas 
swept up by the reverse and forward shocks after the interaction of 
the SN ejecta with the ambient ISM. Thus, the grain size distribu- 
tion calculated by N07 is regarded as that ejected from SNe II to 
the ISM. 

N07 treated two extreme cases for the mixing of elements in 
a SN II: one is the unmixed case in which the original onion-like 
structure of elements is preserved, and the other is the mixed case 
in which the elements are uniformly mixed within the helium core. 
In this paper, we adopt the unmixed case, since the extinction fea- 
tures of carbon and silicon, which are major grain components in 
the unmixed case, are con sistent with observations (Hirash ita et al.1 
l2005l : iKawara et alj 2010). Even if the mixed case is adopted, the 
grain size distribution is similarly biased to large grain sizes, so that 
the behaviour of the extinction curve calculated later is expected to 
be similar. 

As a representative progenitor mass, we adopt 20 Mq, fol- 
lowing N07. The formed grain species are C, Si, Si02, Fe, FeS, 
AI2O3, MgO, MgSi0 3 , and Mg2Si04. According to their calcula- 
tion, small grains with a < 0.02 u.m are trapped in the shocked 
region because the deceleration rate o f a grain by the gas drag 
is inv ersely proportional to its size (e.g. iNozawa. Kozasa. & Habd 
2006). Thus, these small grains are efficiently destroyed by thermal 
sputtering if the ambient hydrogen number density, hh, is larger 
than 0.1 cm -3 . Moreover, the destruction efficiency depends sen- 
sitively on riH- With nu as large as 10 cm -3 , only a few percent 
of grains survive and the grain size distribution is strongly biased 
to large (a > 0.1 tun) radii. It is interesting to point out that 
ISlavin. Jones. & Tielens! fc004l) derived similar grain radii for the 
dynamical decoupling of grains from the interstellar gas, although 
they include magnetic fields in their calculation (see also Section 
|4.4| (. The following results on the effects of grain shattering by tur- 
bulence in WIM are not largely affected as long as the initial size 
distribution is biased to large (o ~ 0.1 [im) grains. 



2.2 Normalization of grain size distribution 

The total grain mass density integrated for all the size range, pdust, 
is normalized to the gas mass density, lAnurriH, to obtain the dust- 
to-gas ratio, T> = pdust/(l-4nHfriH), where mu is the mass of 
hydrogen atom and the factor 1.4 is the correction for the species 
other than hydrogen. Since dust grains are composed of metals, it 
is useful to label the dust abundance in terms of metallicity. We pa- 
rameterize the dust abundance by the oxygen abundance, because 
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Table 1. Characteristic quantities for SN II dust 
production for various ambient densities njj ■ 



n H (cm 3 ) 


m d (Mq) 


■Do 


0.1 


0.35 


2.2 x 10~ 3 


1 


0.14 


8.7 x 10~ 4 


10 


0.033 


2.0 x 10~ 4 



Note, raj is the dust mass ejected per SN II (a 
progenitor mass of 20 Mq is assumed), and T>q 
is the dust-to-gas ratio at the solar metallicity 
(oxygen abundance) Z O 0. 

oxygen is one of the main metals produced by SNe II and oxygen 
emission lines in the optical are often used to estimate the gas-phase 
metal abundance. The oxygen mass produced by a SN II of 20 Mq 
progenitor is mo = 1.58 Mq dUmeda & Nomotdl2002T) . The dust 
mass after the destruction in SNRs (md) is listed for each ambient 
density in Table [T] The solar oxy gen abundance is assumed to be 
Z oQ = 9.7 x 10" 3 in mass ratio flAnders & Grevesselll989h . (The 
oxygen abundance is denoted as Zo and is simply called metallic- 
ity in this paper.) Therefore, in the solar metallicity case, we as- 
sume that the dust-to-gas ratio is T>o = Z Qmd/mo, which is 
also listed in Table [T] The dust-to-gas ratio is assumed to be pro- 
portional to the metallicity: T> — (Zo/Z q)T>q. This is equiva- 
lent to the assumption that both dust and oxygen are predominantly 
supplied from SNe II. 



2.3 Shattering 

Dust grains suffer shattering if the relative velocity between grains 
is larg er than 2.7 and 1. 2 km s for s ilicate and graphite, respec- 
tively i Jones et alii 19961) . As shown bv lYan et al.l J2004l) . these ve- 
locities can be achieved in magnetized and turbulent ISM. In par- 
ticular, grains_are efficiently accelerated in WIM, since damping is 
weak dYan et all 12004). The parameters adopted for MHD turbu- 
lence and the grain velocities obtained are described in Section [2T4l 
The time evolution of grain size distri bution by shattering is calcu- 
lated by adopting the formulation of Jone s et al. Ul994l , ll996l) . We 
briefly review the calculation method. The details are described in 
HY09. 

We solve the shattering equation discretized for the grain size. 
Although nine grain species are treated here (Section l2.lt . the ma- 
terial properties needed for the calculation of shattering are not nec- 
essarily available for all the species. Thus, we divide the grains 
into two groups: one is carbonaceous dust and the other is all the 
other species of dust, and apply the relevant material quantities of 
graphite and silicate, respectively. In fact, as shown later, the mass 
and opacity of the latter group is dominated by Si. The validity 
of this approximation that all the species other than carbonaceous 
dust are treated as silicate (called one-species method) is examined 
in comparison with another extreme approximation (individual- 
species method) in the Appendix [A] Because of the lack of the 
experimental data for Si, we assume that Si (and also the other 'sil- 
icate' species) can be treated as silicate in shattering because of 
similar hardness[f| Since the cratering volume in a grain-grain col- 



lision is approximately proportion al to 1/ P cr , where P CI is the crit- 
ical shock pressure for shattering jjonesetal.il 19961) , the shattering 
time-scale is roughly proportional to P CI . Thus, if P cr is obtained 
for appropriate materials (especially Si) in some future experiment, 
our results can easily be sc aled. The material p roperties of silicate 
and graphite are taken from ljones et al.l j 19961) and summarized in 
HY09. 

The number density of grains whose radii are between a 
and a + da is denoted as n(a) da, where the entire range of a 
is from a m i n to a max . To ensure the conservation of the total 
mass of grains, it is numerically convenient to consider the dis- 
tribution function of grain mass instead of grain size. We denote 
the number density of grains whose masses are between m and 
m + dm as h(m) dm. The two distribution functions are related 
as n(a) da — n(m) dm and m = (47r/3)a 3 p gr , where p gr is the 
grain material density (3.3 and 2.2 g cm -3 for silicate and graphite, 
respectively). 

For numerical calculation, we consider N discrete bins for the 
grain radius. The grain radius in the i-th (i — 1, • ■ • , AO bin is 

— a'^] 1 8, ciQ b ' = o m in, and 



between and a 



where a 



max (i.e. log 5 specifies the width of a logarithmic bin: 
log 8 = (1/N) log(a m ax/a m in)). We represent the grain radius 
and mass in the i-th bin with at = (a^S 1 + a' )/2 and nti = 
(47r/3)af pg r . The boundary of the mass bin is defined as = 
(47r/3)[a^ ] 3 p gr . Giving a m i n , a max , and N, all bins can be set. 
A grain in the i-th bin is called "grain i". We adopt a m i n — 3 x 
10 -4 \jm (3 A) and a max = 3 |J.m to cover the entire grain size 
range in N07, and N = 32. We have confirmed that the results are 
not altered even if N is doubled. 

The mass density of grains contained in the i-th bin, pi, is 
defined as 



Pi = ?nin(mi)(TOj - m i-i) ■ (1) 

Note that pi = pi8i in the expression of IJones et al.l dl994|j l996). 
The time evolution of p; by shattering can be written as 

JV N N 



dpi 



dt 



= —rriip 



j = l k = l 



PkiVki 



if V h i > "shat, 



Qfei = < rmm k 

otherwise, 



(2) 



(3) 



where m^J at (i) is the total mass of the shattered fragments of a 
grain k that enter the i-th bin in the collision between grains k 
and j, <Jki and Vki are, respectively, the grain-grain collision cross 
section and the relative collision speed between grains k and i, and 
fshat i s me velocity threshold for shattering to occur. For the cross 

section, we apply aki = 7i"(ife + at) 2 - 

The grain velocities given bv lYan et alj d2004h are typical ve- 
locity dispersions. In order to take i nto account the dir ectional in- 
formation, we follow the method in I Jones et al. we divide 
each time step into 4 small steps, and we apply Vik = w» + Vk, 
\vi — Vk\, Vi, and Vk in each small step, where Vi and Vk are the 
velocities of grains i and k, respectively (see Section [2~4l l. Note 



3 Amo ng the mate r ials w hose shattering properties are available in Ta- 
ble 1 of Ijones et 

2 Jl996h (i.e. silicate, SiC, ice, iron, and diamond), sil- 
icate is expected to hav e the nearest atomic binding energy. As shown by 
ISerra Dfaz-Cand (2008), who derived the previously unknown shattering 



properties of hydrogenated amorphous carbon, it is not impossible to esti- 
mate relevant quantities, but some guiding quantities from other materials 
is still necessary even in their case. Thus, here we simply adopt the material 
quantities of silicate for Si. 
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that the mass of the shattered fra gment rn*l Ji) depends on Vkj 
as described in IJones etal.l ( fl996h . Briefly, the total fragment mass 
is determined by the shocked mass in the collision and the frag- 
ments are distributed with a grain size distribution oc a~ a{ with 
at — 3.3 unless otherwise stated (see HY09 for the size range of 
the fragments). 

For the shattering duration, several Myr may be appro- 
priate, since it is a typical lifetime of i onizing stars with a 
mass > 20 M (lifetime < 10 Myr) dBressan et al] 1 19931 ; 
llnoue, Hirashita, & Kamavd l2000h . HY09 also indicate that the 
small (a < 0.01 vim) grains depleted by coagulation in dense 
clouds are recovered if shattering in WIM lasts for 3-5 Myr, which 
justifies the necessity of shattering in WIM for mega-years. We also 
examine a longer time-scale, 10 Myr, to investigate a starburst envi- 
ronment where an intense star formation occurs continuously. Such 
a situation may be realized in extragalactic gian t ionized regions 
fairashita & Huntll2004l ; lHunt & H irashita 2009). In summary, we 
examine t = 3, 5, and 10 Myr, where t is the elapsed time of shat- 
tering. 

2.4 Grain velocity 

The grain velocity as a function of grain radius a in the presence 
of interst ellar MHD turbu lence is calculated by the method de- 
scribed in lYan et "all d2004h . They considered the grain acceleration 
by hydrodrag and gyroresonance and calculated the grain velocities 
achieved in various phases of ISM. Among the ISM phases, we fo- 
cus on WIM to investigate the possibility of efficient shattering in 
actively star-forming environments. 

We adopt three cases for the hydrogen number density of WIM 
(n H = 0.1, 1, and 10 cm" 3 ), since N07 applied these densities for 
the ambient medium. For WIM, a density of n u ~ 0.1-1 cm -3 
is usually considered (McKee & Ostriker 1977), but we also ex- 
amine a density as high as riH ~ 10 cm -3 for young H II re- 
gions around massive st ars as observed in starburst environments 
i lHunt & Hirashit j2009h . Embedded starbursts may also have such 
dense ionized regions. We adopt gas temperature T — 8000 K, 
electron number density n c = nu, Alfven speed Va = 20 km 
s -1 and injectio n scale of the turbulence L — 100 pc, following 
lYan et all J2004l) . The effect of the injection scale is minor to that 
of the sound and Alfven velocities. Since both the sound speed and 
the Alfven speed are fixed, the plasma f3 is constant in all cases. 
The g rain charge is assumed to be the same as that in I Yan et"al] 
( 120041) . who calculated it by assuming a typical Galactic condition. 
Since we expect higher interstellar radiation field and higher elec- 
tron density for the starburst environments, the absolute values for 
the grain charge can be larger than those assumed here. For grains 
with a > 0.1 \im, where most of the grain mass is contained in 
our cases, the grain velocity is governed by the gyroresonance. The 
acceleration rate of gyroresonance increases with the grain charge, 
but the accelerati on duration, the gaseou s drag time, decreases with 
the grain charge dYan & Laz arian 2003). As a result, the accelera- 
tion efficiency of gyroresonance is insensitive to the grain charge. 

In Fig. [T] we show the grain velocities. In general, larger 
grains tend to acquire larger velocities because they are coupled 
with larger-scale motions. For small grains, the motion is governed 
by the gaseous drag, w hich has a linear dependence on the grain 
charge dYan et alj|2004h . This is the reason for the complex (non- 
monotonic) behaviour of the grain velocity as a function of a for 
small grains (< 0.1 u.m). We also observe that the grain velocity 
is not very sensitive to tin for large (a > 0.1 u.m) grains, whose 
shattering is important in this paper. 



2.5 Extinction curves 

Extinction curves ha ve been an e ffective tool to examine the 
dust properties (e.g. iMathisI Il990l) . For the calculation of ex- 
tinction curves , we a dopt the same optical constants as those in 
iHirashita et alj ( 120081) for the grain species formed in SNe II (C, Si, 
Si0 2 , Fe, FeS, A1 2 3 , MgO, MgSi0 3 , and Mg 2 Si0 4 ). The grain 
properties of individual species an d the reference s for th e optical 
constants are listed in Table 1 of Hirash ita et al.l J2008b . By us- 
ing those optical constants, we calculate the absorption and scat- 
tering cross sections of homo geneous spherical grains w ith various 
sizes based on the Mie theory (Bohren & Huffma n! 19831) . Then, the 
grain extinction coefficient as a function of wavelength is obtained 
by weighting the cross sections with the grain size distribution. The 
total extinction as a function of wavelength A, denoted as A\, is 
calculated by summing the contribution from all the species. 

As stated in Section 12/31 we divide the grain species into two 
groups in the calculation of shattering: one is carbonaceous dust 
and the other is silicate, which in fact contains all the species other 
than carbonaceous dust. In the calculation of extinction curves, the 
size distribution of the silicate species is redistributed to each com- 
ponent (Si, Si0 2 , Fe, FeS, A1 2 3 , MgO, MgSiOs, and Mg 2 Si0 4 ) 
in proportion to the grain volume (i.e. the total mass of each compo- 
nent divided by its material density) with a fixed shape of the grain 
size distribution. In fact, Si is dominated in the extinction curve, so 
that the uncertainty coming from the above rough treatment does 
not affect our conclusion (Appendix [At. 



3 RESULTS 

3.1 Grain size distribution after shattering 

The grain size distributions after shattering are shown in Figs.[2j]4] 
for riH = 0.1, 1, and 10 cm -3 , respectively. The grain size dis- 
tribution is shown by n(a)/riH- We adopt t = 5 Myr as a typical 
time-scale on which WIM is sustained by the radiation from mas- 
sive stars (Section [2.31 ). Two cases for the metallicity, Zo =0.1 
and 1 Z q, are examined. 

We observe that shattering really affects the grain size distri- 
bution for all the densities. In particular, the abundance of small 
grains with a < 0.1 u.m significantly increases after shattering of a 
small portion of larger grains. If the metallicity is 1 Z q q, a contin- 
uous power-law-like size distribution is realized for a < 0.1 u.m. 
Although the dust abundance is lower for higher nn (Table QJ be- 
cause of more efficient shock destruction in the SN remnant before 
the ejection to the ISM (N07), the grain-grain collision rate is en- 
hanced in higher tih environments. 

The increase of grains with a < 0.1 [im could efficiently af- 
fect the UV and optical extinction curves. This point is quantita- 
tively addressed in Section [3~2l Large grains with a > 0.1 urn 
are marginally affected by shattering; namely, shattering of a small 
fraction of large grains can produce a large number of small grains. 
In the case of HY09, on the other hand, grains with a > 0.1 um 
are more shattered because abunda nt small grains in the MRN 
( iMathis. Rumpl. & Nordsieck 1977]) grain size distribution, which 
they assumed as the initial condition, enhance the grain-grain col- 
lision rate. 

3.2 Extinction curves 

The extinction curve of grains ejected from SNe II tends to be flat 
because small grains are efficiently destroyed in SNRs without es- 
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Figure 1. Grain velocities v calculated from the turbulence model as a function of grain radius a. Two grain species, (a) silicate and (b) carbonaceous dust, 
are shown. The solid, dotted, and dashed lines indicate the velocities with hydrogen number densities of njj = 0.1, 1, and 10 cm -3 , respectively. 




Figure 2. Grain size distributions per hydrogen atom. The solid and dotted lines show the results at t = 5 Myr for metallicities of 1 Z O and 0.1 Z O 0, 
respectively. The hydrogen number density rajj is assumed to be 0.1 cm -3 . The dashed line presents the initial grain size distribution before shattering. Two 
grain species, (a) silicate and (b) carbonaceous dust, are shown. The case with 0. 1 Zqq is multiplied by 10 for the convenience of presentation to offset the 10 
times smaller dust abundance. The arrow is put at a = 0.03 |^m as a rough representative size of the grains contributing to the steepening of the UV extinction 
curve. 



caping into the ISM (Hira shita et alj|2008h . Here we investigate if 
the increase of small grains by shattering effectively steepens the 
extinction curves or not. 

In Fig. [5] we show the time variation of extinction curves for 
rin = 0.1, 1, and 10 cm" 3 with Zo = 1 Z q. W e normalize the 
extinc tion to Ay (i.e. at A = 0.55 |im). As stated by Irlirashita et al.l 
(2008), the initial extinction curve is steeper for lower tih, since 
more small grains survive after the shock destruction in SNRs. We 
also observe that the extinction curve becomes steeper as the grains 



are shattered for a longer time because of the production of small 
grains. Indeed, at t = 5 Myr, A\/Ay at A ~ 0.2 u.m increases 
by more than 20% for hh = 1 cm -3 . The variation of the slope 
by shattering is more pronounced for larger nu since the original 
extinction curve is flatter. 

The extinction curves are dominated by Si and C, which sur- 
vive after the shock destruction in SNRs because of the relatively 
large sizes (N07). Therefore, the steepening of extinction curve is 
mainly due to the production of small Si and C grains by shatter- 
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Figure 4. Same as Fig. 

ing. The contributions from Si, C, and the other species are shown 
in Fig. [6] The 'bump' features at 1/A ~ 4 and 7 |j.m _1 originate 
from the absorption by C and Si, respectively. Such features tend to 
be prominent for smaller grains, since as grains become larger the 
extinction cross sections are mo re determined by the geom etrical 
ones, not by the grain properties teohren & Huffmarj|l9 83). Thus, 
not only the steep slope but also various features in the extinction 
curve become apparent as grains suffer from shattering. 

We also examine the dependence on the dust abundance 




10" 7 10" 6 10" 5 10" 4 
a [cm] 



iut for Tin = 10 cm 3 . 

(metallicity). In Fig.[5jl, we show the evolution of extinction curve 
for Zo =0.1 Z q. The effect of shattering is significantly reduced 
compared with the solar metallicity case. If the grain velocity as a 
function of grain size is fixed, a constant Zot gives the same result. 
Thus, if Zo is ten times lower, 10 times longer time is required for 
the same shattering effect to appear. The scaling with Zot is useful 
if one would like to know the results for other time-scales and/or 
metallicities. 
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Figure 5. Extinction curves normalized to the V band extinction. The solid, dotted, and dot-dashed lines indicate the results at t = 3, 5, and 10 Myr, 
respectively. Panels (a), (b), and (c) present the results for rijj = 0.1, 1, and 10 cm -3 , respectively, with Zq = 1 Zqq. Panel (d) shows the result for 
rin = 1 cm -3 with Zq = 0.1 Zqq. 
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Figure 6. Contributions from Si, C, and the other grain species (dotted, 
dashed, and dot-dashed lines, respectively) for the case of Fig. [5J> (riH = 
1 cm -3 and Zq = 1 Zqq) at t = 10 Myr. The solid line shows the total 
extinction. 



4 DISCUSSION 

4.1 Steepening of UV extinction curve 

In the above section, we have shown that the abundance of small 
(a < 0.1 |J.m) grains increase by shattering of large (a > 0.1 (i.m) 
grains. Consequently, the slope of the UV extinction curve becomes 
steep after shattering. Here, we discuss this issue in terms of the 
grain size distribution. 

The contribution from grains in a logarithmic size range 
[In a, In a + d In o] to the extinction can be written as d/c cxt = 
-Ka 2 Q^(a)n(a)ad\na, where Q\(a) is the extinction cross sec- 
tion normalized to the geometrical cross section (ira 2 ). If the size 
distribution is approximated by a power-law (n oc a~ p ) over a cer- 
tain size range, dK cxt /dlna oc a i ~ v Q\{a). S ince Q\(a) ~ 1 for 
2ira > A and Q\(a) cc a for 2?va <C A (e.g. lBohren & Hu ffman 
1983), we obtain dK cxt /dlna oc a 3 ~ p for 2-na > A and 
dft C xt /d In a oc a 4 ~ p for 2iva -C A. Thus, if p < 3, the largest 
grains have the largest contribution to the extinction. In order for 
small grains to have significant contribution to the extinction, p ^ 3 
should be satisfied. If 3 < p < 4, the largest contribution to the 
extinction comes from the grains with 2na ~ A. In other words, 
the UV (A ~ 0.2 (j.m) extinction curve is steepened significantly 
if grains with a ~ 0.03 [im are produced and p > 3 is satisfied 
around this grain size. 

From Figs. [2]^ we observe that a large number of grains with 
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a ~ 0.03 \xm are produced and the slope around this grain radius 
is p > 3 for the solar metallicity cases. Indeed, the UV slope of ex- 
tinction curve is steepened for the solar metallicity cases as shown 
in Fig. [5] 

In this paper, the shattered fragments are distrib uted with a 
size di stribution with exponent at = 3.3 (Section [2.3| l. |jones et al.l 
from a discussion on the cratering flow, argue that the value 
of at slightly larger than 3 is robust. Even if cif = 2.5 is assumed 
as an extreme case, the difference in the extinction curve is less 
than 10% at A = 0.1 |J.m and smaller at longer wavelengths (see 
the Appendix|B]for details). 



4.2 Grain properties in starburst environments 

From the results above, the presence of small grains in star- 
burst environments is generally predicted, although SNe II tend 
to eject large grains because of the shock destruction in SNRs. 
For example, BCDs (or H II galaxies) in the nearby Universe 
host large ionized region and the age of the current star forma- 
tion episode is a few Myr-20 Myr (e.g. iHirashita & Hunt! [20041 ; 
Take uchi et "all l2005h . These ages are just in the range where 
shattering could modify the grain size distribution and extinc- 
tion curve, although we should take into account the low metal- 
licity in BCDs . Some BCDs show an ex cess of near-infrared 
emission (e.g. iHunt. Vanzi. & Thuarj 1200 lh . which can be at- 
tributed by the emission from transiently hea t ed very small grains 



The time-scale on which shattering in SN shocks effectively 
destroy s large grains can basically estimated by a similar way to 
iMcKeel Jl989h . A single SN can sweep M 8W ~ 10 4 M of gas 
(i.e. .M sw Us/2 ~ Z?sn with a shock velocity v B ~ 100 km s" 1 



1979 



ISellgrenl fl984l ; iDraine & Anderson! 
2005) have carried out a comprehensive 



i Aannestad & Kenvon 
Il985h . iGalliano et al.l 
analysis of the SEDs of dust and stars in some dwarf galaxies 
(dwarf irregular galaxies and BCDs), and have shown that the grain 
size is biased to small grains with ~ a few nm. Since their sample 
galaxies have metallicities larger than 1/10 Zq, shattering in WIM 
can work as a production source of nm-sized grains on time-scales 
of a few Myr and thus can be considered as an origin of small grains 
in these galaxies. 

It is natural to expect that a similar condition (i.e. turbulence 
in WIM sustained more than a few Myr) is generally realized in 
starburst galaxies. Although it is hard to compare the extinction 
curve with the observed wavelength dependence of the dust atten- 
uation beca use of the effects of radiative transfer dCalzettill200ll ; 
Ilnouell2005l) . shattering may be crucial to reproduce the reddening 
in starburst galaxies. Therefore, shattering should be considered as 
a source of small grains, which contribute to the reddening. Or dust 
produced by AG B stars in the underly ing old population (older than 
several x 10 s vn lValiante et alj2009h could contribute to the steep- 
ening if they produce small grains; however, there are some ob- 
servational indicatio ns that dust grains produced in AGB sta rs are 
large (a ~ 0.1 (i.m) JGroenewege n 1997; Gauger et al. 1999). 

Efficie nt shattering also o ccurs in the ISM by the passage of 
SN shocks. Ijones et alj dl996h show that a large fraction of large 
grains with a > 0.1 |J.m is redistributed into smaller grains by 
a single passage of shock with a velocity of ~ 100 km s _1 . In 
their calculation, large grains take longer time before they are dy- 
namic ally coupled with gas and are subject to more collisions with 
dust. Ijones et aL Ul996l) consider the MRN distribution as the ini- 
tial grain size distribution, which enhances the shattering efficiency 
compared with our case, because of the enhanced collision with the 
abundant small grains. Below, we estimate the time- sca le on which 
shatter ing in SN shocks destroys large grains based on I Jones et al.l 
( 1996), although the time-scale obtained might be an underestimate 
for the grains produced by SNe U, because of the enhanced colli- 
sion rate in the MRN distribution. 



and energy given to gas by a SN iJsi- 



10 01 erg). Then, the 



gas mass swept by SN shocks with v s > 100 km s~ per unit 
time can be estimated as Af sw 7, where 7 is the SN rate. Thus, the 
time-scale on which the entire gas mass M s is affected by shat- 
tering by SN shocks is estimated as r sw ~ M g / (M sw 7). Since 
7/V ~ 10 _2 Mq 1 for a Salpeter initial mass function JSalpeterl 
Il955h (ip is the star formation rate), the above time-scale is esti- 
mated as t sw ~ 10 -2 Mg/tp. This estimate indicates that the shat- 
tering time-scale by SN shocks is about 0.01 times the gas con- 
sumption time-scale by star formation. In starburst environm ents, 
M g /V> ~ 10 8 -10 9 yr may be reasonable jYoung et all 19861) . and 
shattering in SN shocks occurs in 1-10 Myr, which is comparable 
to the time-scale investigated in this paper. Therefore, both shatter- 
ing in turbulence and that in SN shocks can affect the grain size 
distribution. A detailed calculation of shattering in SN shocks of 
grains produced by SNe II is required before we judge which of 
these two shattering mechanisms is dominated. 

It might be also useful to discuss our results in terms of the ex- 
tinction curves of the Large and Small Magellanic Cloud (LMC and 
SMC), both o f which have develo ped H II regions such as 30 Do- 
radus. Indeed. [Bernard et alj fc008l) indicate that the 70 \xm excess 
around 30 Doradus can be explained by an enhancement of the 
abund ance of very sma ll grains possibly by the destru ction of large 
rains . iBot et alj [2004) find this excess in the SMC. IParadis et al.l 
20091) show that the very small grain abundance is really enhanced 
around 30 Doradus by using an SED model of dust emission. How- 
ever, the extinction curves in these galaxies are much steeper than 
our results (A\/A v ~ 2.9 and 3.2 a t A ~ 0.2 ^.m for the LMC 
and the SMC, respectively; |Pei|[ 1992 b- Since those galaxies have 
less intense star formation than BCDs, it is hard to extract the star- 
bursting components where shattering of large grains should be 
working as investigated in this paper. The steep extinction curves 
of the LMC and the SMC indicate that we should consider not 
only the dust production/shattering in star-forming regions but also 
some other mechanisms which act as efficient production sources 
of small grains. For example, shattering in warm neutral medium 
works on a time-scale of 100 Myr (HY09). ISM phase exchange, 
which occurs on a time-scale of 50-100 Myr, also affects th e evo- 
lution of grain size distribution jO'Dormell & Mathisll 19971) . Such 
longer-time-scale mechanisms could also have affected the extinc- 
tion curves (grain size distributions) of those galaxies. The current 
paper, which focuses on a short-time-scale (< 10 Myr) grain pro- 
cessing, is a starting point to include other physical processes in 
future work. 



4.3 Comparison with high-z data 

At 2 > 5, it is usually assumed that the main production source 
of dust is SNe II whose progenitors have short lifetimes, since 
the cosmic age is too young for low mass stars to evolve (but 
sec IValiante et aI1l2009l) . Thus, the extinction curves at such high 
z are often used to te s t the theory of dust production in SNe II 
dMaiolino et aT]|2004bl ; IHirashita et aT]|2005h . As a representative 
case of observed high-z extinction curve, we discuss the restframe 
UV extinction curve o f SDSS J 1048+4637 (z = 6.2) obtained by 
iMaiolino etafl d2004bl) . 

In Fig. [7] we show the UV part of the extinction curves calcu- 
lated by our models in comparison with the observed UV extinction 
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Figure 7. The same extinction curves as those shown in Fig.[3p are plotted 
only in the UV range (Zq = 1 Zqq and njj = 1 cm -3 ). The shaded area 
shows the observed ex tinction curve for SDSS J1048+4637 (z = 6.2) by 
Maio lino et aT] (2004b), including the uncertainty. 

curve of SDSS J1048+4637. The extinction curves are normalized 
to the value at A = 0.3 \im. We show the result for riH = 1 cm -3 , 
but the following disc ussions hold qualitativ ely also for other den- 
sities. As discussed in lHirashita et al.l j2008h . the initial extinction 
curve before shattering is too flat to explain the UV rise in the 
observed extinction curve because small grains are selectively de- 
stroyed in SNRs. However, after shattering, the extinction curve 
approaches the observed curve because of the production of small 
grains. After 10 Myr of shattering, the observed extinction curve 
is reproduced. Not only the slope but also the bump feature at 
1/A ~ 4 vim , which becomes prominent after shattering (Sec- 
tion ^. 2t , may account for the behaviour of the observed extinction 
curve around 1/A ~ 3.5—4 \im~ . 

In summary, if the metallicity is nearly solar and the age of 
the current episode of starburst is larger than 5 Myr, we should 
take the effect of shattering in turbulence into account in comparing 
the observed extinction curve with the theoretical one even at z > 
5. Since quasars tend to be found in evolved stellar s ystem whose 
metal licity could be nearly solar (or more than solar; ljuaerz et al.l 
2009), the UV rise of the extinction curve may be caused by the 
production of small grains by shattering. The dependence of the 
extinction curve on age and metallicity may also be responsible for 
the variation of UV sl ope of the quasar spectra in the sample of 
iMaiolino et ai]b004al) . 



netic fields in shocks produce complexity in the kinematics of large 
(> 0.1 irm) grains. Thus, it may be important to trace the grain tra- 
jectory around the reverse and forward shocks. The quantification 
of all these effects of magnetic fields is left for future work. 

Nevertheless, the importance of shattering by turbulence for 
small-grain production in starburst galaxies should be an important 
issue even if we consider the the effect of magnetic field in the fu- 
ture, because it is still true that the shock destruction in SNRs sup- 
presses the injection of small grains into the ISM. It should also be 
kept in mind that at the smallest size ranges (a few A), the treatment 
of grains as bulk solid may not be a good approximation. Since 
such tiny grains do not affect the UV-optical extinction curve as 
discussed in Section |4~T1 the results on the extinction curves are not 
affected. Mid-infrared spectra of dust emission are more suitable to 
constrain the abundance of such small grains (e.g. lMathiJl990h . 



5 CONCLUSION 

We have theoretically investigated the effect of shattering in turbu- 
lent WIM on the grain size distributio n by using the framework for 
shattering bv ljones et all ( I1994I.1 19961) and th e calculation of inter- 
stellar MHD turbulence by Yan et al. (2004). We have focused on 
systems in which dust is predominantly produced by SNe II. Al- 
though SNe II tend to eject large (a > 0.1 irm) grains because of 
the shock destruction in SNRs (N07), shattering in WIM supplies 
small grains on a time-scale of several Myr in the solar-metallicity 
(i.e. Galactic dust-to-gas ratio) case. Consequently, the extinction 
curve is steepened and the features such as the carbon bump around 
1/A ~ 4 fim" 1 and the Si bump around 1/A ~ 7 [im -1 become 
apparent if the metallicity is solar and the duration of shattering is 
longer than ~ 5 Myr. Therefore, when we treat a system in which 
the metallicity is solar and the star formation age is > 5 Myr, we 
should take into account the effect of shattering in interstellar turbu- 
lence. In particular, the extinction curves of high-z quasars, whose 
metallicity is typically (above) solar, may be affected by shattering, 
and the UV rise of the extinction curve as well as the bump fea- 
ture at 1/A ~ 3.5—4 p.m can be attributed to the small grains 
produced by shattering. If the metallicity is < 1/10 solar, the ex- 
tinction curve does not vary significantly on a time-scale of < 10 
Myr because the frequency of grain-grain collision is reduced in 
proportion to the grain abundance. Thus, the steepening mecha- 
nism of extinction curve discussed in this paper is valid for systems 
whose metallicities are significantly larger than 1/10 solar. We con- 
clude that shattering in WIM is generally of potential importance 
in starburst galaxies as a production mechanism of small grains. 



4.4 Remarks on grain physics 

Before concluding this paper, we mention some physical pro- 
cesses to be considered in the future. In the calculation of 
the shock destruction of grains in SNRs by N07, the ef- 
fects of grain electrical charge an d the effect s of magnetic 



1996) 



2007) 



fields are ignored. As shown in Ijones et all dl994l 

and more recently by iGuillet. Pineau Pes Forets. & Jonesj i 

and IGuillet. Jones. & P ineau Des Foretsl J2009I) . the dynamics of 
charged grains is critically modified by magnetic fields. The gy- 
ration around the magnetic fields tend to strengthen the coupling 
between gas and dust, and this effect could suppress the ejection of 
large grains into the ISM. Thus, not only small grains but also large 
grains with a > 0.1 yim could b e subject to significant processing 
in the shock. Isiavin et "all j2004l) show that the presence of mag- 
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APPENDIX A: TEST FOR THE ' ONE-SPECIES ' METHOD 

As stated in Section l2~3l all the grain species other than carbona- 
ceous grains are treated as a single species, called formally 'sili- 
cate' in calculating the grain size distribution. This approximation 
is called 'one-species' method, and it is exact if all the grain species 
have the same shape of grain size distribution. We expect that the 
one-species method gives a reasonable answer since Si is domi- 
nated among the 'silicate' category. Although the 'silicate' species 
other than Si (we call these species non-Si grains) have minor con- 
tributions in grain mass, some of them have a significant contribu- 
tion to the number of small-sized grains, which affect the UV slope 
of the extinction curve. Here we test the validity of the one-species 
method in comparison with the 'individual-species method' as ex- 
plained below. 

The 'individual-species' method adopts the grain size distribu- 
tion of individual species and the evolution of grain size distribution 
is separately calculated for individual species (note that the grain 
size distribution summed over all the species other than carbona- 
ceous grains is adopted for 'silicate' in the one-species method). 
In calculating the evolution of grain size distribution of a certain 
species, the total mass density of the species relative to the gas 
density is assumed to be the total dust-to-gas ratio (but the grain 
size distribution after shattering is normalized again to recover the 
correct mass ratio of each species). This treatment maximizes the 
production of small grains for non-Si species, which have smaller 
sizes than Si, but minimizes the production of small Si grains. Thus, 
this method is suitable to examine the maximum possible contribu- 
tion from non-Si small grains to the UV extinction curve. 

In Fig. IAU we compare the grain size distributions predicted 
by the one-species and individual-species methods for nu = 0.1 
and 1 cm~ 3 at 5 Myr. For nu = 10 cm -3 , the difference between 
the two methods is negligible because non-Si grains contribute little 
to the total grain abundance. From the figure, we observe that the 
difference is relatively large in the case of nu —0.1 cm -3 . This is 
because the fraction of non-Si grains is larger for nu =0.1 cm~ 3 
than for nu = 1 cm -3 . 
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Figure Al. Size distributions for (a) njr = 0.1 cm -3 and (b) rijr = 1 cm -3 for the grains other than carbon (i.e. 'silicate'). The solid and dotted lines show 
the results with the individual-species method and with the one-species method, respectively. The dashed line presents the contribution from Si to the solid 
line. The metallicity and the age are assumed to be 1 Zqq and 5 Myr, respectively. 
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Figure A2. Extinction curves normalized to the V band extinction for the grain size distributions presented in Fig. lAll fthe grain size distributions of carbona- 
ceous grains are the same as those in Figs.|2]and[3j. Panels (a) and (b) present the cases with njj =0.1 and 1 cm - 3 , respectively. The solid and dotted lines 
show the results of the individual-species and one-species methods, respectively. The dashed line represents the contributions from various species (C, Si, and 
the others as labelled in the figures) for the individual-species method. 



In Fig. IA2I we show the extinction curves calculated by 
the two methods. We observe that the extinction curves of the 
individual-species method tend to be steeper than those of the one- 
species method. As can be seen in the figure, the steeper slope 
comes from the contribution from the non-Si grains indicated by 
'others'. In the individual-species method, the size distribution of 
each non-Si species, which has a larger fraction of small grains 
than that of Si, is calculated separately, so that the production 
of small non-Si grains is enhanced. We note that the 'real' grain 
size distribution would lie between the results of the two methods. 
This means the approximate treatment adopted in the text (i.e. one- 
species method) is justified for nrr > 1 cm" 3 . For hh < 0.1 cm -3 , 



because the contribution from non-Si species is significant, the er- 
ror of the one-species method is at most ~ 10% at A ~ 0.2 \im, 
and ~ 40% at A ~ 0.1 \im. In order to overcome this uncertainty, 
we should develop a different scheme that could treat the collisions 
between multiple species (in our case, 9 species), which the current 
scheme cannot treat in a reasonable computational time. 
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Figure Bl. Grain size distributions for njj = 1 cm -3 with a metallicity of 1 Zqq. The solid and dotted lines show the results at t = 5 Myr for Of = 3.3 
and 2.5, respectively. The dashed line presents the initial grain size distribution before shattering. Two grain species, (a) silicate and (b) carbonaceous dust, are 
shown. The arrow is put at a = 0.03 nm as a rough representative size of the grains contributing to the steepening of the UV extinction curve. 



APPENDIX B: FRAGMENT SIZE DISTRIBUTION WITH 
A SHALLOWER SLOPE 

The size distribution of shattered fragments is assumed to be a 
power law with an exponent of — cvf , As we discuss in the text, the 
steepening of extinction curve becomes prominent if the power-law 
exponent (p) of the grain size distribution around a ~ 0.03 \xm is 
steeper than ~ 3 (Section [4. lt . ljones et al.l (fl996) have shown that 
the size distribution after shattering is not sensitive to at . They also 
argue that at slightly larger than 3 is robust against the change of 
the cratering flow parameters in shattering (ctf = 3.3 is adopted in 
the text). Nevertheless it would be interesting to examine if p > 3 
is realized even if we assume at < 3. 

Here we examine the smallest exponent adopted in ljones et al.l 

dl996h . Qf = 2.5 as an extreme case. The ambient hydrogen num- 
ber density is fixed to nu = 1 cm -3 . In Fig. lBll we show the re- 
sult at t = 5 Myr. As expected, the effect of Qf is more prominent 
for smaller grains, since shattering with large ctf can supply small 
grains more efficiently. However, we observe that the difference be- 
tween Qf = 2.5 and 3.3 is sma ll around a ~ 0.03 \im, confirming 
the result o f lJones et alj Jl996l) . The small difference comes from 
the fixed shattered mass in a collision; that is, the distribution of 
grain fragments as a function of size has a minor effect compared 
with the total mass of shattered fragments (shattering efficiency). 

The extinction curves are shown in Fig. IB2I We observe that 
the difference between the two curves with a — 2.5 and 3.3 is 
negligibly small at A ~ 0.3 [im and is less than 10% even at 
A ~ 0.1 \im. The small difference is the natural consequence of 
the small variation of grain size distribution at a > 0.03 \im. 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 
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Figure B2. Extinction curves normalized to the V band extinction for the 
grain size distributions presented in Fig. lBll The solid and dotted lines show 
the results for ot{ = 3.3 and 2.5, respectively. The dashed line represents 
the initial extinction curve before shattering. 



